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Sub Decoherence Time Generation and Detection of Orbital Entanglement 

F. Brange, O. Malkoc, and P. Samuelsson 
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Recent experiments have demonstrated sub decoherence time control of individual single-electron 
orbital qubits. Here we propose a quantum dot based scheme for generation and detection of 
pairs of orbitally entangled electrons on a timescale much shorter than the decoherence time. The 
electrons are entangled, via two-particle interference, and transferred to the detectors during a single 
cotunneling event, making the scheme insensitive to charge noise. For sufficiently long detector dot 
lifetimes, cross-correlation detection of the dot charges can be performed with real-time counting 
techniques, opening up for an unambiguous short-time Bell inequality test of orbital entanglement. 

PACS numbers: 73.63.Kv, 03.65.Ud, 03.67.Bg, 73.50.Td 


The concept of quantum entanglement has ever since 
its inception attracted much attention. Initially ques¬ 
tioned because of its nonlocal properties, violating local 
realism 0 , 0 , entanglement has over the past decades 
emerged as an indi^ensable resource for quantum infor¬ 
mation processing [3| . Spurred by pr^osals for electronic 
spin-based quantum computing [J, , spin qubit exper¬ 

iments 0 , 0 and demonstrations of long spin decoher¬ 
ence times [ 8 | , large efforts have been devoted to investi¬ 
gations of spin entanglement in nanostructures. Recent 
experimental progress comprises entanglement of single¬ 
electron 0 and two-electron [l^ spin qubits and splitting 
[lll - [l^ of spin-singlet Cooper pairs in hybrid supercon¬ 
ducting systems. 

In contrast to spin, entanglement between electronic 


orbital degrees of freedom [ 1 ^ [ig , such as charge states 
in quantum dots 17[ or edge channels in quantum Hall 
systems [isl - li^ , has received limited attention. In partic¬ 
ular, orbital entanglement has not been demonstrated ex¬ 
perimentally. The key reason is arguably that superposi¬ 
tions of orbital states are sensitive to charge noise, result¬ 
ing in short decoherence times, of the order of nanosec¬ 
onds 2ll-24|. This has led to the widespread view that, 
despite all-electrical quantum state control and read-out, 
electronic orbital degrees of freedom cannot be harnessed 
for quantum information processing. Very recently this 
view was contested by demonstrations of fast, coherent 
operations of single-electron orbital qubits on the pi¬ 
cosecond timescale [25l - E^ . several orders of magnitude 
shorter than the decoherence time. These experiments 
motivate renewed efforts on orbital-based quantum infor¬ 
mation processing and call for novel schemes to generate 
and detect orbital entanglement on timescales well below 
the decoherence time. 

Here we propose such an entan glem ent scheme, based 
on coherent electron cotunneling (28l| in a quantum dot 
system, see Fig. [T] During the cotunneling event, of 
the order of picoseconds 29, [s^, the electrons are en¬ 
tangled via two-particle interference Bin and simul¬ 
taneously transferred to the detectors, fully preserving 
coherence 3^ ^. We show, based on the full transfer 
statistics 34]35|, that the entanglement can conveniently 


be detected by violating a Bell inequality (BI) formu¬ 
lated in terms of low-frequency current cross-correlators 
liBB. Moreover, for long enough detector dot life¬ 


times, measurements of coincident electrons in the detec¬ 
tor dots can be performed with real-time charge counting 
techniques [37|. This opens up for an unambiguous BI- 
test of orbitally entangled electrons in nanosystems based 
on short-time measurements. 



FIG. 1: (a) Schematics of the setup, with entangler {E) and 
detector {A and B) subsystems, quantum dots (circles), elec¬ 
tronic leads (hatched rectangles), and key tnnneling couplings 
(thick lines) shown, (b) A typical configuration of E and A, B 
dot energies ei, £ 2 , Co, ep with a = A±, P = BE and inter-dot 
charging energies C/ 12 , Uap shown, fulfilling the two-particle 
resonance condition ei + €2 + C /12 = Ca + ep + Uap. (c) The 
three tunneling processes in the pair transfer cycle: (I) En¬ 
tangler dots 1, 2 populated with rates Fi, r 2 . (II) Coherent 
cotunneling of electrons from entangler to detector dots (here 
a = A— and /3 = B+), with amplitude tpa 2 i- (HI) Electrons 
in a, P emitted to detector leads, with rates F^, Up- 


The combined entangler-detector system, shown in 
Fig. [U consists of two entangler {E) dots, 1 and 2, and 
four detector dots, a = A± at A and /3 = BE at B. 
Each dot 7 = 1,2, a, pi has one active, spin-degenerate 
level at energy and is tunnel-coupled to an electronic 
lead, with a rate F-^. The six dots are further coupled to 
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each other (nearest neighbors) with tunneling amplitudes 
tai,ta 2 ,tf 3 i,ti 32 ,ti 2 ,tA+A- andts+B-- The Coulomb in¬ 
teraction energy between different dots, 7 ^ 7 ', is Ujj>. 
Due to strong on-site repulsion, double occupation of the 
dots is prevented. Consequently, the spin degree of free¬ 
dom only leads to a renormalization of tunneling rates 
and is hereafter neglected. 

As illustrated in Fig. [U the dot level energies are 
tuned to suppress single-particle tunneling between en- 
tangler and detector dots (as well as between the two en- 
tangler dots). In addition, to make cotunneling the domi¬ 
nating entangler-detector transport mechanism, the ener¬ 
gies are chosen to optimize the conditions + ~ 

Cl + £2 + U 12 for resonant two-particle tunneling between 
entangler dots 1,2 and detector dots a, /3. The amplitude 
for the cotunneling is denoted Moreover, cotunnel¬ 

ing between the entangler dots and two dots at the same 
detector, A or B, is tuned off resonance. 

The dot-lead tunneling rates obey ra,r^ ::!> ri,r 2 
t,Sa 2 i/^- Hence, the cotunneling is much slower than 
the dot-lead tunneling (sequential) and back-tunneling 
from the detector to the entangler dots is suppressed. 
Moreover, in this regime there are at most two particles 
in the dot-system at the same time. Furthermore, taking 
hr a ;§> tA+A- and hFfj ^ ts+B-, tunneling between 
the detector dots can be neglected. The two entangler 
leads are kept at a finite bias while the four detector 
dots are grounded. Throughout the paper we consider 
the high-bias regime; the energy difference between dot 
system states differing by one electron is well inside the 
bias window and the lead temperature can effectively be 
put to zero. 

We first give a physically compelling picture of the 
ideal working of the entangler-detector system, a detailed 
analysis follows below. In Fig. [U the tunneling processes 
in the pair transfer cycle, from entangler to detector 
leads, are shown. (I) Starting from an empty system, 
electrons tunnel sequentially into the two entangler dots 
1 and 2. (II) The electrons in the entangler dots cotunnel 
resonantly to the detector dots, one electron to a and one 
to (3. This process is coherent and key to the entangler- 
detector scheme in the following ways: (1) The electrons 
can cotunnel in two ways, from I to a and 2 to /3 or 1 
to P and 2 to a. This gives rise to orbital entanglement 
via two-particle interference as illustrated in Fig. [2j The 
quantum state emitted from the entangler can be written 

|^)ent = C 2 i| 2 )^| 1 )b — Ci 2 | 1 )a| 2 )_B ( 1 ) 

with | 1 )^ denoting an electron emitted from dot 1 to¬ 
wards detector A etc. and C 12 , C 21 constants depend¬ 
ing on the system properties. We note that for |ci 2 | = 
|c 2 i| = |'0)ent is maximally entangled. (2) By tun¬ 

ing, with electrostatic gates, the dot-entangler tunneling 
amplitudes tia etc., the emitted state |' 0 )ent can be lo¬ 
cally rotated in the | 1 ),| 2 ) basis during the transfer to 
the detector dots, see Fig. O 


(III) At the detector dots, the two particles tunnel se¬ 
quentially out to the leads. As will be described below, 
depending on the values of F^j, F^ the entanglement can 
be detected via violation of a HI formulated in terms of 
either low-frequency current cross-correlations or short- 
time measurements of joint detection probabilities. 



FIG. 2: Schematics of the entanglement creation and state ro¬ 
tation during the cotunneling, emphasizing the orbital states 
|1)a|2)_b (|2)a|1)s) for electrons emitted towards A (B). The 
state |i/’)ent in Eq. [T]is a coherent two-particle superposition, 
describing the entangler dot electrons being emitted 1 —>■ A, 
2 ^ B (filled dots, state |1)a|2)s) or 2 — >■ A, 1 ^ i? (empty 
dots, |2)a|1)s). The tunability of the tunneling amplitudes 
tai etc., here illustrated as beam splitters (shaded yellow) pa¬ 
rameterized by angles 6 a , 63, allows for a rotation in orbital 
{|1), |2) |-|-), |—)}-space of A and B. 


To substantiate the qualitative picture above we first 
recall that uncertainty principle arguments give a cotun¬ 
neling time ^ h/AE, where AE is determined by the 
energy of the classically forbidden state, virtually popu¬ 
lated during the tunneling [ 2 ^. Here we take the cotun¬ 
neling time, typically of order of picoseconds in quantum 
dot systems [29| , much smaller than the decoherence time 
l/F,^, of the order of nanoseconds [U [ 2 ^ . The cotun¬ 
neling process is thus coherent, as recently demonstrated 
for double-dot systems [s^ , and can be treated within an 
effective Hamiltonian approach. 

To this end, we write the total Hamiltonian as E[ = 
i?o + H. The term i?o describes the single-particle dot 
levels, the Coulomb repulsion between the dots, the elec¬ 
tronic leads and the dot-lead tunneling. The perturba¬ 
tion V describes the dot-dot tunneling. The cotunneling 
dynamics is obtained from a Schrieffer-Wolff transforma¬ 
tion [^, where V is eliminated to leading order, see [s^ . 
Besides cotunneling processes, the transformation yields 
renormalization terms which are absorbed into the dot- 
level energies, the Coulomb interaction strengths and the 
dot-lead couplings in iJo- Furthermore, neglecting terms 
describing processes much slower than competing sequen¬ 
tial dot-lead tunneling or resonant two-particle cotunnel¬ 
ing, we obtain the effective cotunneling Hamiltonian 


HeS = Ho+^ (ti3a2idldld2di + H.c.^ 

0.(5 


( 2 ) 


with djj {dp) the creation (annihilation) operator of an 
electron in dot p. Here the two-particle tunneling ampli¬ 
tudes tpa 2 i are related to the single-particle amplitudes 
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through 


, _ ^/31^ck2 ^al^/32 /q\ 

with the energies AEpa = AE^pla ■(->■ /3] depending on 
e-y and Uj^/, see [s^. In particular, at two-particle reso¬ 
nance, €a. E €/3 E UafS = Cl -|- £2 + Ui 2 = E^, 1/ AEp^ = 
{Er — [^2 + £/3 + U 213 ]) ^ + {Er — [ci + €a + Ula]) ^■ 

Importantly, the terms in Eq. (jS]) correspond to two 
different two-particle paths, I —>■ /3, 2 —>■ a and 1 ^ a, 
2 ^ P (see Fig. [I, that take the two electrons be¬ 
tween the entangler and detector dots, with tunneling 
times Ti/AEPa and Ti/AEap^ respectively. Hence, the to¬ 
tal cotunneling amplitude tpa 2 i is a coherent superposi¬ 
tion of the individual two-particle tunneling amplitudes. 
The minus sign between the terms is a consequence of 
fermionic exchange. Note that the amplitude for both 
electrons to tunnel to the same detector side, A or B, is 
negligible since the process is off-resonant (due to level 
detuning and Coulomb interaction Ua+a-,Ub+b-)- 

Starting from the state with both entangler dots oc¬ 
cupied, time evolution governed by iJeff gives that the 
two-particle state tunneling out into the detector dots 
is \i>)dot OC with |a) = |-|-)^,|-)a and 

1 ^) = |+)b, |-)s- To arrive at the result in Eq. o , we 
focus on the experimentally relevant regime where the 
specific dot (± in a, /3) dependent part of the Coulomb 
interactions between detector dots at A and B, Uap, 
as well as between entangler and detector dots, Uia 
etc., is negligibly small compared to the relevant energy 
level differences in AEap, AEpa- At two-particle reso¬ 
nance, this regime implies Ca ^ ca, £/3 ^ and thus 
AEap = AEab, AEpa = AEba, independent of a,/3 
(see Fig. [T|). In addition, we consider the single-particle 
tunneling amplitudes to be tuned to fulfill the relation 
\tA+iltA-i\ = \tA- 2 /tA+ 2 \- Under these conditions, with 
Eq. (|3]), one can write the state in the detector dots 

I*)-- = S. = ( “"/j- ) (4) 

with Si the scattering matrix for the effective beam split¬ 
ter at detector i = A, B, see Fig. [5] Here we have in¬ 
troduced the parameterization of the tunneling ampli¬ 
tudes tA+i = IaiCOsOa, tA -1 = tAishlOA ctc., with 
|tAiP = |tyi+iP + etc. Moreover, we can directly 

identify ci 2 = tAitB 2 / AEab and C 21 = tA 2 tBi/ AEba 
in Eq. (HD. Importantly, Eq. o verifies the claim that 
the cotunneling from the entangler to the detector dots 
can be described as an emission of an entangled state, 
|■^/')ent, locally rotated in orbital space before arriving to 
the detector dots. 

Detection of the entanglement of |'(/')ent can be per¬ 
formed by transport measurements. In the high-bias 
regime considered, the full transport statistics of the 
entangler-detector system can be described within the 


framework of a Markovian quantum master equation 
^ = C-yy[p] for the reduced density operator p, with the 
Liouvillian superoperator 


Ex[p\ — ^Ws 1 p] + p] 


+ (1 - f^)Vy,^[d^,p]) + ^V0[dld^,p] 


( 5 ) 


where Hs is the Hamiltonian of the dot system, is the 
lead Fermi function with fj = l (0) for 7 = 1,2 (a,/I), 
V^lL^p] = e^'^LpL^ — i {L' 1 L,p} is the dissipator, and 
X-y is a counting field for lead 7 . The last term in Eq. ([5]) 
describes dephasing with a rate independent of 7 . 

The cumulant generating function Fy, is obtained as 
the eigenvalue of fulfilling the condition lim^^o ^ 
0. To leading order in the tunneling amplitudes tpa 2 i we 
find 

Fy, = J2 (e*A<«+x/3-xi-X2) _ (g) 

Q(/3 


where the transfer rates Pap are given by 


Pap = 


|t/3a21 P (Tct + r^ d-E,^) 

X + E/J + E,^)^ -f ef2a/3 


( 7 ) 


withei2a/3 = eo,-l-e/3-l-C/a/3-[ei-l-e2-l-C/i2] the energy away 
from two-particle resonance. The generating function 
in Eq. ([SD, a sum of terms — gives 

a physically clear picture of the full transport statistics: 
There are four independent, elementary events of electron 
pair transfer from the entangler leads 1, 2 to the detector 
leads a and /3 with rates Pap/^h. 

To formulate a Bl-test in terms of transport quanti¬ 
ties, we consider either the dephasing-broadened regime, 
where E,^ ^ Eq,,E^, or the situation where Eq,,E ;3 are 
independent of the specific dot (± in a, P). In both cases 
Pap fx |t/ 3 Q; 2 iP, with a proportionality constant indepen¬ 
dent of a, p. This allows us to identify, directly from 
Eq. ([6]), Pap/YlapEap as the probability that a pair of 
electrons emitted from E is jointly detected, with one 
electron in lead a and one in p. Based on these joint de¬ 
tection mobabilities, we can directly formulate a CHSH 
BI [illii as 


S = \Eab - Ea'b + Eab' + Ea'b' \ < 2 , ( 8 ) 


with the correlation functions Eab = 
Pa+b+-Pa+b--Pa-b++Pa-b- A,A',B,B' 

Fa+b++^a+b--\-^a-b+-^^a-b- ’ ’ ’ 

denote different detector settings, i.e., here different 
single-particle tunneling amplitudes. Importantly, the 
rates Pap can be obtained from the current cross- 

correlation Sap = lx=o = Pa/3, allowing for 

an experimental test of entanglement via a violation of 
the BI. We note that for the state \'tjj)ent, the maximal 
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Bell parameter is S'max = 2 \/l + sin^ 9 [l^, where 
6 = 2 arctan(ci 2 /c 2 i), i.e., all entangled states can in 
principle be detected by a BI violation. 

Beyond being insensitive to decoherence, does our 
setup offer any additional advantages, not present in ex¬ 
isting proposals? The answer is yes. By tuning the de¬ 
tector tunneling barriers such that the rates Tq,, be¬ 
come small, typically in the sub-MHz range, the tunnel¬ 
ing on and off the detector dots can be monitored in real 
time, via, e.g., time-dependent electrical currents flowing 
through quantum point contacts capacitively coupled to 
the dots (not shown in Fig. [IJ [44-l46|. In particular, by 
simultaneously detecting the c harg e transfers at dots a 
and /3, as demonstrated in Ref. [37| . one can identify the 
coincidence probabilities for individual, emitted pairs ar¬ 
riving at the detector dots or leads. This opens up for 
an entanglement test based on a violation of a BI formu¬ 
lated in terms of short-time joint detection probabilities, 
in direct analogy to quantum optics. 

Theoretically, the short-time properties of the corre¬ 
lated charge transfer can be described by the electronic 
analogue [i^ of Glauber’s second degree of coherence [i^ 
in quantum optics, 


= 


{{JqH(t)J;3)) -F {{Jpn{T)Ja)) 
2((Jc))((T/3)) 


(9) 


where ((A)) is the stationary expectation value of oper¬ 
ator A, H(r) = 6^°'^ is the master equation propagator 
and is the current superoperator at 

lead j = a, p. We note that the Liouvillian is given 
by Eq. ([S]) and that, by definition, e{{Ja)) = la- To 
leading order in tpa 2 i, we obtain 




TaT^ 


■ Tfl 2I„L 


a-l/3 


■(e 


-r„r -^g-rpr) 


( 10 ) 


Several important conclusions can be drawn from this 
result. First, » 1 for all times t < l/Tcl/T/j 

(terms in g\p{T) of order unity are neglected), a typi¬ 
cal signature for pair transport [i^. Second, the decay 
timescales of the correlations, I/Tq, 1/Tp are set by the 
tunneling times out of the dots a and /3. Third, g^^^ (r) 
by construction describes the correlations between elec¬ 
trons arriving at the leads a, /3, i.e., tunneling out of the 
corresponding dots a, /3. To obtain the probabilities that 
two electrons arrive coincidentally at the detector dots a 
and j3, we have to integrate g^pi''') times, i.e., 

accounting for all possible emissions from the dots to the 
leads. This gives the quantity (not normalized with cur¬ 
rents) 


21alp 


J drg^^liT) = e^Pafi 


( 11 ) 


which thus gives the same probabilities Pap / 

as in the long-time limit, Eqs. m and 0, producing a 

consistent theoretical picture of the charge transfer. 


Interestingly, in contrast to charge detection via long¬ 
time current correlations, real-time detection of coin¬ 
cident probabilities has typically non-ideal efficiencies, 
Va.gp < 1- The efficiency at A is given by r]a = 

(and similarly at B), where Aw is the bandwith of the 
detector. Based on existing experiments, near-unity effi¬ 
ciency of the detection requires TcT^ well below 1 MHz 

ii-lll. 


In conclusion, we have proposed a quantum dot based 
entangler-detector system that generates and detects or- 
bitally entangled electrons on a timescale much shorter 
than the decoherence time. The main idea is to use co¬ 
tunneling for both the entangling and state manipulation 
during the transfer to the detectors. Recent experimental 
demonstrations of highly controllable multiple quantum 
dots systems [IB-il together with existing short-time 
charge detection techniques 44-4^ make the realization 
of our scheme within experimental reach. 
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Here we derive the effective cotunneling Hamiltonian 
in Eq. (2) in the main text. We start by considering 
the isolated dot-system described by the Hamiltonian 
Hs = Hq + V, where Hq = ^ , 

7 77' 

with single-particle dot-levels and the Coulomb re¬ 
pulsion between different dots (7 ^ 7 '). Here 

n.y = d\d^ is the number operator. The term V — 

S M^77' dt^dj/ + t*Ydyd.y) describes the tunneling, with 
77' 

amplitudes between nearest neighbor dots. We set 
the tunneling amplitudes between detector dots effec¬ 
tively to zero as tunneling out to the leads will dominate 
over back-tunneling to other dots. For the remaining 
single-particle dot-dot tunneling couplings we assume off 
resonance conditions, allowing us to eliminate V to first 
order in Hs by means of a generalized Schrieffer-Wolff 
transformation Hs —>■ e~^Hse^ = Heg [l|. 

The generator S of the transformation is obtained from 
the condition [l| 

[S,Ho] = -V (1) 

To find an explicit expression of S, we first introduce an 
operator 


A. 


77' 


-Sbrv- 


+ E(1^7S'-17ys') 




6'eB 


seB sgb 


that [Ajji,Ho] = 0 , since both operators only consist of 
number operators. Importantly, we have 


^A^ryf d\.d^', Hq^ — Ay^'^dtyd^', Hq^ 


= A 


77 


Cj A 'y u^§t^Ti§f 


d^ d^/ 



n n 

seB 5eB 


dt^d-yl = —dl^d-y! 


where we note that the sum on the second line is over all 
dots except 7 and 7 '. We then get 


[S,H^] = 

1 


yy/ 

^ ^ 2 (^77^ [^77^ ^7^7^ ’ "^0] ^yy' [^7yC?y/C^7; 

^ ^ “ ^tyy/(i^(iy' “t“ ^,-yyfdy/d'J^ - V 


which shows that S in Eq. ([5)) fulfills Eq. (Ill- 

Next, we compute the effective Hamiltonian to leading 
order in V from [i| 

HeB = Ho + ^[S,V] (3) 


for 7 7 ^ 7 b To clarify the notation, we point out that 
A^^i contains number operators acting only on the four 
dots (of the in total six) which are not 7 , 7 '. Out of these 
four dots, B denotes a configuration of n = 0—4 occupied 
dots and B denotes the remaining 4 —n unoccupied dots, 
i.e., the complement of B. The sum in A^^i runs over all 
16 possible configurations B. Moreover, each term in the 

sum contains a denominator e-,, —e.y'-|- ^ {U^S' — 

S'^B 

the energy difference associated with single-particle tun¬ 
neling between 7 ' and 7 in presence of the occupation 
configuration B of the four other dots. 

In terms of A^y ; f generator S can then be expressed 
as 

^ = E ^A7'(77'4rf7' - tl^Y^X^d^) (2) 

77 ' 

To prove this, we note that A^y is Hermitian and ful¬ 
fills the property A~^y — —A^y- Furthermore, we note 


To evaluate the commutator of S and V, we first consider 
the commutator of a single term of S, Ass'itss'dgds’ — 
iggidgids), and a single term of V, tjjid'jydji + t*^,dl^,dj. 

Four different terms are produced by the commutator, 
dependent on the relation between 7 and <5 and between 
7 ' and S' (note that 7 7 ^ 7 ' and 6 7 ^ S'). Below we 
evaluate and discuss them one by one, making use of the 
assumption that there can be maximally two electrons in 
the dot system at the same time. 

(i) <5 = 7,<5' = y 

We have the commutator 


rj [-^77' 1^77' l.yry/d.yf dn^) , d.yd.y' t / d f dy) 


= |b 


yY Ayyl {Uy — Uyl) (4) 


This term only consists of number operators and can 
therefore be eliminated by means of a renormalization 
of the single-particle energy levels Cy and the interaction 
strengths Uyy'. 
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(ii) 5 = 7 , 5' 7 ^ 7 ' and (in) <5 7 ^ 7 , 5' = 7 ' 

For (5 = 7 , (5' 7 ^ 7 ', we have the commutator 

(i'yS't>y§fd^fd'y'jj d)ydj/ “t“ try^fd^td^^ 

- \n^l ^ ^ 7 ] (tyS' ^yy' dry/ dg/ -\- try^'tyy/ dg/ d^> ) (S) 

where A^g/ [n^/ n^] denotes A-yS ', with n^/ replaced by 

Ury. This term describes tunneling of a single particle, 
first from 5' to 7 and then from 7 to 7 ', or vice versa. A 
similar term is obtained when 5 ^ 5' = . 

However, in the considered parameter regime (speci¬ 
fied in the main text), these tunneling processes are not 
important for the dynamics of the dot system, as there 
is always a competing process which dominates. To be 
explicit, there are three different processes. First, if an 
electron occupies one of the detector dots, the most likely 
process is tunneling to the lead and not tunneling to an¬ 
other dot, since the dot-lead couplings are much stronger 
than the dot-dot couplings. Second, if one of the entan- 
gler dots is occupied, the most likely process is that the 
second entangler dot becomes occupied. Third, when 
both entangler dots are occupied, the dominating pro¬ 
cess is cotunneling, see (iv), since this process is close to 
resonance in contrast to single-particle tunneling. The 
single-particle tunneling processes described by Eq. m 


are thus neglected. 

(iv) d ^ 7 ,d' 7 ^ 7 ' 

We have the commutator 

A( _!_!_) 

2 \e5 -I- Ujs — £6’ — Urys/ es -I- Ury/s — es’ — Uj/s' J 

X ^sS't /yyfdgd^^ d/y' d^^ A ^SSA^y^/dg/dry/dryds 

A ^SSAyy' d^f d)ydryf ds -t- tgg/tryry/d^dry/drydg/'^ (6) 

These terms describe two-particle cotunneling processes, 
e.g., the term oc dgdtydj/ds/ describes tunneling from dots 
7 ', S' to 7 , S. Importantly, only the terms describing co¬ 
tunneling between the two entangler dots 1 , 2 and two 
detector dots a, jd are relevant for the dot system dy¬ 
namics, all other cotunneling processes have competing 
single-particle tunneling processes which dominate, sim¬ 
ilar to (ii) and (hi). 

The total two-particle tunneling amplitude t/ 3 a 2 i be¬ 
tween the entangler and detector dots is given by sum¬ 
ming up all terms in S and V, yielding the relation 
tj 3 a 2 i = (^) main text, where 


I/AA^q, — „ ( 


1 


1 


1 


2 Vei -I- Ul2 — e/3 — 1^2/3 £2 + 1^2/3 — Ca — Uaji (-2 + Ul2 — ta — Ula Cl -|- Ula — 6/3 — UaP 


(7) 




1/1 1 1 1 
2 Ve 2 -f Ui 2 — €/3 — Uip ei -I- Uip — Cq, — Uap ei -I- U 12 — ia — U 2 a £2 -|- U 2 a ~ £,S ~ Uafi 

I- 


) ( 8 ) 


Denoting with Ei, Ef and Em the initial, final and 
intermediate energies, respectively, of the system during 
the two-particle cotunneling, 1/Ai5^„, lAEap can each 
be written as \ J2m( Ei-Em ~ Em-Ef where m runs 
over the two possible intermediate states. We note that 
on two-particle resonance, Ef = Ca + ep + Uap = £i + £2 + 
U 12 = Ei, this simplifies to e -e ’ ™ accordance 
with standard second-order perturbation theory Q. 

Using the two-particle tunneling amplitude in Eqs. © 
and ([5|), the effective Hamiltonian can be written as 

iJeff = Ho + (ti3a2ldldl^d2di + (9) 

q/3 

which is our Eq. (2) in the main text. 

We finally note that a full treatment, including the 
terms describing the leads and the dot-lead couplings 
in Hq (not presented here), gives an additional level¬ 


broadening term in the denominator of A^^ry/, propor¬ 
tional to the coupling strength F^, F^ between a lead 
and its dot. Since these coupling strengths are much 
smaller than the differences between the single-particle 
energy levels of the dots, this level-broadening effect is 
negligible. Furthermore, cotunneling processes involving 
leads are suppressed in the regime considered, where the 
dot-lead couplings are weak compared to the energy dif¬ 
ferences between single-particle energy levels of the dots. 
Hence, to leading order, the effect of the leads can be 
neglected in the Hamiltonian. 
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